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Abstract 

We deal with a problem of finding maximum of a function from the Holder class 
on a quantum computer. We show matching lower and upper bounds on the 
complexity of this problem. We prove upper bounds by constructing an algorithm 
that uses the algorithm for finding maximum of a discrete sequence. To prove 
lower bounds we use results for finding logical OR of sequence of bits. We show 
that quantum computation yields a quadratic speed-up over deterministic and 
randomized algorithms. 

1 Introduction 

Quantum algorithms yield a speed-up over deterministic and Monte Carlo algorithms 
for many problems. Many papers deal with quantum solution of discrete problems, 
starting from the work of Shor followed by database search algorithm of Grover 
0. Other discrete problems were also studied, such as discrete summation, computa- 
tion of the mean, median and fcth-smallest element |3*].[1].|5].[7].|12|. 
There is also a progress in studying the quantum complexity of numerical problems. The 
first paper dealing with a continuous problem was that of Novak |14j . who considered 
integration of a function from the Holder class. The problem of function approximation 
on quantum computer was studied by Heinrich |E1,[S]- Also path integration |16j and 
differential equations ^0])^B on a c[ uan t um computer were investigated. 
In this paper, we deal with a problem of finding maximum of a function from the Holder 
class on a quantum computer. The complexity of this problem in deterministic and 
randomized settings on a classical computer is well known 13 . 

We present matching upper and lower complexity bounds in the quantum setting. We 
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show that quantum computations yield a quadratic speed-up compared to determin- 
istic and randomized algorithms over entire range of class parameters. Upper bounds 
are shown by constructing a suitable algorithm, which uses the optimal algorithm for 
finding maximum of a discrete sequence. To prove the lower bound we use the result 
of Nayak and Wu 12 for finding logical OR of sequence of bits. 

In the next section necessary definitions are presented. Existing results for searching 
maximum of a sequence are shown in Section El The main result of this paper is 
contained in Theorem ^ in Section 

2 Quantum setting 

In this section we briefly describe the model of computation. More details about quan- 
tum computing can be found in Jl]. Suppose that we have a numerical problem given 
by a solution operator S : F — > G, where F is a subset of a linear function space and G 
is a normed space. We are interested in computing an approximation of S(f) for / 6 F 
on a quantum computer. This is done by an algorithm A. The algorithm can access 
the input element / only by a quantum query. An output of the algorithm A for a 
given / is a random variable A(f, u). For a detailed discussion of quantum algorithms 
and the quantum query operator, the reader is referred to |14j . 

We know recall what is meant by the error of an algorithm. Let < 6 < 1. The local 
error of algorithm A on input element / is defined by 

e quant (5)A/j0) = ^ . P Q S (f) - A(J,U,)\\ > e} < 9}. 

The number 1 — 9 is thus the lower bound on success probability. For e > the bound 
e quant A, f, 9) < e holds iff algorithm A computes S(f) with error at most e and 
probability at least 1 — 9. The global error in the class F is defined as 

e qUant (S, Aj F; 0) = gup e quant ( ^ ^ ^ ^ 

For 9 = 1/4, we denote the error for / by e quant (S', A, f) and the error in the class F 
by e quant (5, A,F). 

The cost of the algorithm on the input element /, cost (A, f) is defined as a number of 
accesses to an oracle. In classical settings it is a number of function values or derivative 
values that is used to compute maximum of the function. In the quantum setting, by 
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an oracle we mean the quantum query operator. The global cost of algorithm is defined 

as 

cost(j4, F) = sup cost {A, /). 
feF 

For e > 0, the e-complexity is defined as the minimal cost of an algorithm that produces 
an e-approximation: 

compf ant (5,F) =min{cost(^,F)| e quant (S, A, F) < e}. 

In the next section some known results about complexity of finding maximum of a 
sequence of numbers in the quantum and classical settings are presented. 

3 Searching maximum of a sequence 

We recall results on a discrete maximization (minimization) problem. Consider the 
following problem: given a sequence X = (xq,x±, . . . , x n _i) of real numbers in [0, 1], 
find the number Xi = max(min){a;,- : j = 0, 1, . . . , n — 1}. The cost of an algorithm is 
defined as a number of accesses to the oracle, which returns the input number x^, i = 
0, 1, . . . , n — 1. Another possibility is to count a number of comparisons. In this model 
the oracle returns the logical value of a comparison Xi < Xj, where i,j 6 {0,1,... , n— 1}. 
Clearly, the complexity of this problem on a classical computer is 

comp worst (n) = 9(n) , comp rand (n) = 9(n), (1) 

where by comp(n) we mean the minimal cost of an algorithm computing the maximum 
(minimum) value from the sequence of n numbers in suitable setting. Better results 
can be obtained on a quantum computer. In 1996 C. Diirr and P. Hoyer in [S] pre- 
sented comparison quantum algorithm for finding the minimum. This algorithm finds 
minimum value from the list of n items with probability greater than 1/2 and its run- 
ning time is 0(y/n). They based their algorithm on quantum exponential searching 
algorithm 2 , which is a generalization of Grover's search algorithm introduced in [gj. 
This result establishes the upper bound on the complexity of the problem of finding 
maximum of a discrete sequence on a quantum computer. 

Lower bounds on this problem were established by A. Nayak and F. Wu in ^2]. They 
examine the more general problem: for X = (xq, ,xi, . . . ,x n -i) G [0, l] n and A > 
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compute A-approximate fcth-smallest element, i.e., a number X{ that is a j'th-smallest 
element of X for some integer j 6 (k — A, k + A). 

If A = 1 (or less) this problem reduces to the problem of finding /cth-smallest element 
exactly. For k = n — 1 fcth-smallest element is the maximum value from the sequence. 
In ^2], a quantum algorithm has been presented with the cost 0{N\og(N) log log (N)), 



where N = yn/A + \Jk{n — k)/A. The algorithm is inspired by the minimum finding 
algorithm of Diirr and H0yer 0, and uses exponential search algorithm of Boyer et al. 
PJ. It finds A-approximation of £;th-smailest element for any k E {0, . . . , n — 1} and 
A > 1 with probability at leat 2/3. 

Nayak and Wu in established essentially matching lower bounds for this problem. 
To derive the bounds they used polynomial method introduced by R. Beals et al. in 

HI- 

These results show that the complexity of searching maximum of n elements on a 
quantum computer is of order 

comp quant (n) = 8(Vn). (2) 

The comparison of this result with Q shows that quantum computers make a quadratic 
speed-up over classical computers for this problem. 

4 Searching maximum of a function 

We consider the problem of finding the maximum of a function from the Holder class 

F^ = {/:[0,l] d ^R| f€C r , H/ll <1, 

D^f{x)-D^f{y)\<\\x-yf Vx,yG[0,l] d , V£>W}, 

where run through the set of all partial derivatives of order r, r G No, < p < 1 
and || • || = || • ||oo- We want to find a number 

M(/)= max f(t) 
te[o,i] d 

up to some given precision e > 0, for any function / from class F^' p with probability 
not less than 3/4. 

The complexity of this problem on a classical computer in the deterministic and ran- 
domized settings is presented in [13]. We shall recall these results for a further com- 
parison. 
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In the deterministic worst-case setting the local error of algorithm A on input function 
/ £ Fj P is defined by 

e ™*( M ,AJ) = \M(f)-A{f)\ 

and the global error by 

e wmst {M,A,F r d ' p ) = sup e worst (M,yl,/). 

In the randomized setting, algorithm A = (A(u>)) u ^ is a random variable on some 
probabilistic space (Q,B,m). The local error of this algorithm is defined by 

e rand (M,A/)= / \M(f)-A(u)(f)\dm(<j), 
Jn 

and the global error by 

e TanA {M,A,F r /) = sup e™ nA {M,A,f). 

The cost of an algorithm in the deterministic and randomized settings is meant as a 
number of function values accessed by an algorithm. In the randomized setting, points 
where / is evaluated can be chosen randomly. 

It is shown in ^3] (PP- 34 and 59) that the complexity of function maximization in the 
Holder class in both worst-case and randomized settings is given by 

( d \ / d \ 

) and comp^ d (M,F^) = 9 ' [ [ V 

We now pass to the quantum setting. In the following result we prove that a significant 
improvement is achieved on a quantum computer. 

Theorem 1 Let e > 0. The quantum e-complexity of computing maximum of a func- 
tion from Holder class F,' p is 

compf ant (M,Fj' p ) = 6 (V^T 
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Proof: 

First we prove the upper bound. We divide each edge of the cube [0, l] d into n subin- 
tervals of equal length. We get N = n d cubes K l , i = 1, . . . , N. Let t l denote the center 
of cube K l . On every cube K l we use Taylor's expansion of /. For t 6 K l we have 



where 



and 



Let 



w 



f(t)=w i (t) + R r (t,t i ), 



i (*) = Erf/ (fc) (* i )(*- ti ) fc 5 

fc=o K - 



R r (t, f) = jf 1 (f^(et + (i - e)f) - )) (t - ey !LJLl de. 



m i(f) = max w % {t). 
We consider the algorithm A* of the form 

A *(f) = . ™ mi(/), 

8=1,..., N 

where fhi(f) is an approximation of rrii(f) computed by some classical algorithm. We 
assume that 

\Mf)-Mf)\<£l Vi = l,-..,N, 

for some E\ > independent of i and /. To compute fhi(f) on a classical computer 
we do not need any new evaluations of / or its partial derivatives, so that information 
cost does not increase. (Of course, it is still not an easy task to compute fhi(f) and 
it increases combinatory cost of the algorithm.) The maximum of the discrete set of 
numbers fhi(f), . . . ,fhjsf(f) we compute on a quantum computer, with probability not 
less than |. This is done by the optimal algorithm described in Section 13 
We now estimate the error of the algorithm defined above 

e^ at (M,A*J) = \M(f) - A*(f)\ = | max /(i) - . max rfn(f)\ 

te[o,i] d 8=1,..., N 

= I max max f (t) — max fhAf)\ 
1 i=l,...,N tzK* v ; i=l,...,N u 71 

< max I max f(t) — fhi{f)\ 
i=l,...,N t€K l 



< max ( | max/(t) - m»(/)| + |mj(/) - m»(/)| I 

i=l,...,N V t&K 1 J 
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= max I max f (t) — max w % (t) I + ei 
i=i,...,jv ieK< teic* 

< max max|/(t) - wUt)] + e x . (3) 

i=l,...,N t€K l 

For i £ i^ 1 we have 

\f(t)-w\t)\ = \R r (t,f)\ 

C + (1 - 0)t< ) - )) (t - * i ) r (1 ( ~_ g) 1 ^ 1 dg 

< sup \f^{et + {i-e)t)-^\f)\ \\t-tf f 1 ^f^— do 

< h sup \\et + (i - e)t - tY \\t - 1* ir 

0e[o,i] 

= H \\t - f\\ r+p < H ( ■ 

The constant H depends on d and r but not on n. From this, and inequality (j3J) we 
get the error bound 

e quant (M,^*,/) <h(-Y P + ei. 



We now choose £\ = (l/n) r+p . Then for some constant G independent of n we have 

e quant (M, A*,f) < G {^j +P • (4) 
We examine the cost of algorithm A*. 

To compute fhi(f) we need to know the value of / and the values of all its partial 
derivatives of order up to r at point t l . So the cost of computing fhi(f) is 

'd + k-l\ (d + r)l 



d\r\ 



which is independent on n. 

Due to (J2J the cost of computing the maximum of numbers fhi(f), • • • , fhi(f) on quan- 
tum computer is 0(^/N) = 0(VnM) accesses to the numbers rhi(f). Thus, the total 
cost is 

cost (A*,/) < Cn d ' 2 (5) 
for some constant C. Due to (jlj), to obtain e qua,nt (M , A* , f) < e it suffices to take 



K e 2 ( r +p) function and derivative values, where K is a constant. This completes the 
proof of the upper bound. 
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We now prove the lower bound. Assume that A is any algorithm that computes 
M(f) = max tg [ 01 ]d f(t) for any / £ F^ p up to the error e with probability not less 
than |. We denote the cost of this algorithm by c(e). 



For si > 0, the class F r ^ p contains n = ye l r+p J functions /i, • • • ,/ n with disjoint 
supports such that max tg [ 01 ]d fi(t) = e\ (see jT3], p. 35). 

Let E\ = Ae. Let X = (x±, . . . , x n ) be any sequence such that Xi G {0, 1} Vi = 1, . . . , n. 
Then the function 

n 

/ei := ^ ] S-i/i 
j=l 

belongs to the class F^' p . Thus, algorithm A applied to f £1 computes max tG j 01 ]d f £l (f) 
up to the error e = with the cost c(s) = c(x)- That is 

| max / ei (t)-A(/ ei )|<5 (6) 

with probability not less than |, and cost c(e) = c(^-). 
From the definition of f £l we see that 

, I £ i if maxj=i n xi = 1 

max / £1 (t) = < . (7) 

te[0,l] d if max i= i ... n Xi = 



If A(f £l ) > |, then due to © 



1 1 



So, in this case, due to Q, max^m yd = £i and max^i,...^ Xj = 1 with probability 

at least |. Similarly, with probability at least |, if A(f si ) < \e\, then max tg [ 01 ]d f £l (t) = 
and maxj=i „ X{ = 0. 

Based on algorithm A, we now define an algorithm A, which finds the maximum of a 
sequence X = (x%, . . . , x n ). This algorithm is constructed as follows: 

if |ei < A(f £l ) < |ei , then we put A(X) = 1, and 
if ~\ei < A(f £l ) < \e x , then we put A(X) = 0. 

In the other cases we put A(X) = 0. With probability at least |, we have that 
maxj=i n Xi = and A(X) = 0, or maXi=i j ... n x^ = 1 and A{X) = 1. Hence, the 
algorithm A computes the ixLcixiniuni of ti numbers xi, . . . , x n ^ such that xi S {0, 1} 
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(logical OR of the input bits), with probability not less than |. From (see proof of 
Theorem 1.5 in the case of Xj 6 {0, 1}) we know that the cost of such an algorithm is 
Q(-y/n). The cost of the algorithm A is not less than the cost of the algorithm A. Since 



n 




we have that 



c(e) = (Vn) = n ( Q) 
This completes the proof of the theorem. 



d 

1 \ 2(r+p) 



Comparing this to the classical deterministic or random complexity of this problem, 
which is G , we see that a quantum computer makes a quadratic speed-up 

over classical settings. This is achieved over the entire range of r, p and d. For the 
integration problem, a quadratic speed-up over the randomized setting holds only for 
(r + p) jd small. 
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